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Abstract 

We apply Coupled Cluster Method to a strongly correlated lattice Hamil- 
tonian and we extend the Coupled Cluster linear response method to the cal- 
culation of electronic spectra. We do so by finding an approximation to a re- 
solvent operator which describes the spectral response of the Coupled Cluster 
solution to excitation operators. In our Spectral Coupled Cluster Method the 
ground and excited states appear as resonances in the spectra and the resol- 
vent can be iteratively improved in selected spectral regions. We apply our 
method to a Mn02 plane model which corresponds to previous experimental 
works. 

1 Introduction 

The numerical methods for solid state physics span a wide range of tech- 
niques which aim to provide approximate solutions to the problem of many 
body interactions in correlated systems; the exact solution being unknown. 
Between these techniques some of the most notables are: the DMF theory, 
which uses a self-energy correction term obtained from an Anderson impu- 
rity many-body solution, the GW equation, which calculates the self-energy 
with neglection of vertex corrections, the quantum monte-carlo method and 
the Coupled-Cluster Method (CCM). The Coupled-Cluster method has been 
conceived in the 5Qs by Fritz Coester and Hermann Kummel for nuclear 
physics, and has since progressively gained other domains |1 1. In quantum- 
chemistry, in particular, CCM is widely regarded as the most reliable choice 
when high accuracy is needed|2|. Concerning the lattice models of strongly 
interacting electrons CCM has recently been applied to spin lattices[3| and 
to Hubbard model|4|. 

Although the CCM was initially formulated as a ground state approxima- 
tion, the recent developpment of CCM time dependent linear response |5 1 has 
stretched the CCM applicability to excited states. In particular Crawford and 
Ruud have calculated vibrational eigenstates contributions to Raman optical 
activity|6| while Govind et al. have calculated excitonic states in potassium 
bromide Q. 
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In this paper we extend the coupled cluster linear response method to the 
calculation of electronic excitation spectra. To do so we represent an initial 
wave-function as the product of the probe operator times the CCM solution, 
and we develop an original solution method of the resolvent equation for the 
CCM ansatz. 

This paper is organised in the following way. In section |2] we detail the 
equations. In section |3] we describe the model of a Mn02 plane derived 
from previous absorption and scattering x-ray spectroscopy studies, on which 
we test our method. We discuss the results of our Spectral Coupled Cluster 
Method in section |4] There we also validate the method by comparing it 
to the exact solution that we can obtain when we restrict the Hilbert space 
dimension to such an extent that the exact diagonahsation is possible. 



2 Method 

In the Coupled-Cluster methodll] one searches an approximated solution to 
the eigen-problem 

(1) 

where H is the Hamiltonian in second quantization and is formed by a sum 
of products of one-particle creation and annihilation operators. The one- 
particle operators change, between 1 and 0, the occupation integer number of 
the one-particle orbitals contained in the model. 

The solution is represented, given a reference state |4>o), by the exponen- 
tial ansatz 

1^) = lOo) ~ lOo) (2) 

where S is the ideal exact solution and Sn is a sum, truncated to terms, of 
products of electron-hole pair excitations : 

Sn^Y^ ti Symm I Y\ c^^^cl^ , \ (3) 
1=1 [<.-=i ' ' J 

In this formula is the number of degrees of freedom of the ansatz. The 
larger is this number the more accurate is the representation. The r, 's are free 
coefficients that must be obtained from the CCM equations below. Each term 
in the sum is the product of a set of electron(hole)-creation operators . Each 
term is determined by a choice of indexes a, a, ,t), with the greek(latin) 
letter a(a) ranging over empty(occupied) orbitals. One can cosiders the ref- 
erence state as the vacuum state and that each term / in the sum Sn creates, 
from vacuum, an excited state which is populated by n, particles (holes and/or 
electrons). The Symm operator makes the ansatz symmetric for the Hamilto- 
nian symmetry subgroup which transforms, up to a factor, the reference state 
|Oo) into itself. 

In the Coupled Cluster method, the rationale for the exponential ansatz 
resides in its size extensivity property. This means that for a system com- 
posed of two non correlated parts, A and B, the coupled cluster ansatz oper- 
ator can be factorized as the product of two operators e^^^^ = e^^e^^ . This 
simple factorisability relation has deep consequences| 1 1, whose one of the 
most important is that, in a system with periodic translational symmetry, the 
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calculation complexity for a given accuracy does not depends on the system 
size. 

The CCM equations are obtained substituting j*?), in the eigen-equation, 
with its ansatz and by multiplying at the left with e^^'^ , the inverse of the 
ansatz operator . One obtains for the eigenvalue 



while the free parameters are obtained setting the eigen-equation residue to 
zero in the space of excited states which enter the Sn sum : 



The Coupled Cluster Method expands these equations by means of the 
Hausdorff expansion formula which for two arbitrary operators A and B states 
that: 



e-^Be^ =B+[B,A] + 1/2[[B,A],A] + ..l/n![[B,A]...],A] + ... (6) 



The numerical applicability of CCM relies on the fact that, when A is substi- 
tuted with and H replaces B, only the first five terms in the series, can be 
non zero. This can be demonstrated considering that Sn is formed by creation 
operators only, and that the interactions contained H are composed by prod- 
ucts of up to four single-particle operators for the Coulomb interaction. For 
each term of the expansion every Sn entering in the commutators must have 
at least one one-particle creation operator contracted with one annihilation 
operator of H, for the term not to be identically zero. 

Equation |5] gives polynomial equations by which we can determine 
the N unknowns f,. These equations have order up to the fourth in the f, 
variables because this is the maximum order in Sn for the non-zero terms of 
the Hausdorff expansion. 

The number of solutions of a system of polynomial equations explodes 
exponentially with the number of equations and it is not possible, except for 
small systems, to explore systematically the whole solutions space. To solve 
the equations, instead, we use the Newton's method to follow the solution, 
increasing iteratively the number of free parameters and using as a starting 
point the for parameters the — 1 ones found at the previous iteration plus 
a random choice for the A^''' one. 

The accuracy of the CCM solution increases with A^. At each iteration 
the new (A^+ 1)'^' term is constructed, in equation [s] by assigning its order 
fi{N+i) ^nd by choosing the concerned electron and hole orbitals which are 
expressed by the sets of indexes OC(Ar+i) and af^N+\).h with k ranging from 1 
to ni^N+i) ■ We denote the ensemble of all possible choices with the symbol 



(4) 




(5) 




(7) 



where the possible choices satisfy the condition 




(8) 
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The simplest choice consists in choosing the ^ which gives the largest residue 
in equation |8] 

Once we have obtained the CCM ground state and its ground energy E, 
we are interested in the transition probability for a time dependent pertur- 
bation exp{i(j3ot)D, where D is an arbitrary product of operators. The 
transition rate is given by the Fermi golden rule which states that the proba- 
bility for the absorption of an energy quantum/zcOo =fi<X) — E, with/zco being 
the final state energy, is proportional to: 

. , _ {<Po\e^'D^H-&-iy)-^De^\<i>o} 

Pd((0,y)=3w^ — ^ : — - (9) 

(4>o|e'^'e'^|4>o> 

where y is a small line width. In order to calculate the above expression we 
have to solve two problems : find an approximate solution R for the resolvent 
equation : 

{H-(o-iy)\R>=De^\^Q) (10) 

and calculate the scalar product. 

We represent an approximated solution for the resolvent, introducing the 
approximating operator Ro.a.y and the following ansatz which is similar to 
the ansatz for S with the difference that it contains both annihilation and 
creation operators and that, in order to accede to the whole spectra, no sym- 
metrization is done : : 

\R>= RD.co.yDe^l^o) (11) 

RD,(i>.y= ''DM.y + ^f=l''D.m,yl\k'=l^ji.k ^^^"^ 

in this expression r' are free parameters and we have introduced the nota- 
tion c to represent in a compact way both creation and annihilation operators. 
The definition of the c operator is, naming by N^rbs the total number of rep- 
resented orbitals (occupied and empty) : 



^ J-Norbs J ^ i^orbs + 1,2* Norbs] 



(13) 



We build our spectral CCM equations (SCCM equations) by multiplying 
at the left with e^^'^, and by setting the residue to zero : 

1= {<Po\D^e-^{H-Oi-iy)RD,^.yDe^\<pQ) (14) 

0= {^o\D'' (nlliC;„)%-^(//-CO-/Y)/?D.a,.YZ)e^|4>o) V/ G [1,A^'-](15) 

Note that the validity of the above equations relies on the fact that D, being 
a product of operators, commutes with S. These equations are expanded 
by the Hausdorff expansion formula substituting, in equation [6] A with Sn 
and B with (// — co — i'Y)RD.<>i.y- The Hausdorff expansion contains, also in 
this case, a finite number of non-zero terms because each term of the the 



resolvent operator , in equation 12 contains a finite number of annihilation 
operators while, as discussed above, each term of H contains a maximum of 
four annihilation operators. 
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The expansion gives a set of linear equations for the r parameters. The 
resolvent equation accuracy is improved by systematically increasing A^'', se- 
lecting, at each iteration, the set of numbers 



{(" 



iV'+l ' V'+l.O' V'+lji^,.^l J J ^^^^ 

corresponding to the largest residue in the SCCM equations. When we calcu- 
late the residue we fix co = co, at the center of the spectral region of interest. 
Over the spectral region of interest the r parameters are given by a linear 
algebra operations of the kind r = (Mi)^^ {M2 + C0M3) where the M's are 
matrices obtained from SCCM expansion. Once we know the R operator we 
can calculate the spectra with the following equation : 

This expression can be expanded using the Wick's theorem and the linked- 
cluster theorem as akeady done by Sourav et al.[8J. Contracting in all possi- 
ble ways the operators contained in D^Ro^a.yD with themselves one obtains 
sum of products of Green's function of different orders. To simplify this 
we use the simplest approximation which consists in setting to zero all the 
connected Green's function excepted the one-particle Green's function : 



GUiJi)^ ' " , i'\ (18) 



We expand this equation for G using the Wick's and the linked-cluster theo- 
rems. We obtain a hierarchical set of equations involving Green's functions 
of arbitrary order. This expansion needs to be truncated choosing a closure 
relation. This closure relation is already provided by the choice that we have 
made setting to zero all the connected Green's functions except the two points 
one. To obtain the Dyson equation for the Green's function we proceed in the 
following way. Each time we contract a c operator with one of the terms con- 
tained in S, on the right, or with 5^ on the left, a new vertex is obtained from 
which a number of new lines, equal to the term order minus one, are coming 
out. We consider all the combinatorial ways of contracting all these lines 
with themselves, except one branch which propagates further the Green's 
function. This is the analogous of the Hartree-Fock approximation where 
two of the four legs of each Coulomb vertex are contracted with each other. 
The Dyson equation is solved iteratively. The final result for the spectral 



function of equation 17 depends linearly on the parameters r, functions of 



CO, and contain products of Green's functions ( defined by equation 18 1. The 
spectral resonances positions depend on the r parameters, which are found 
by the SCCM equation and whose behavior accounts for many-body corre- 
lations. The resonances intensities, instead, depend on our Hartree-Fock-like 
truncation which still accounts many-body interactions but in the mean-field 
approximation. 
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3 Model 



In previous studies on manganites we applied exact diagonalisation, and 
Lanczos method, to the study of resonant X-ray scattering ||9l at L2X3 edges 
and A'p fluorescence flOl. The spectroscopy data were modeled with a small 
planar cluster, described in second quantization. The model consisted of the 
central Mn atom open shells orbitals, plus some selected orbitals localised on 
the first neighbouring shell of oxygens atoms and Mn atoms. These studies 
revealed a pronounced O 2p character of the doped charge carriers, and the 
non-local nature of the forces governing the charge redistribution phenom- 
ena which are very important in these systems. The accounting of few extra 
orbitals from neighbouring shells, beside the resonating atom, is crucial in 
describing these phenomena but one rapidly ecounters the limit of the ex- 
ponential growth of the Hilbert space dimension, when trying to extend the 
size of the cluster. To calculate the ground state and spectra of larger sys- 
tems, while still keeping a good description of the many-body correlations, 
we have developed the methods described in this paper 

We will compare SCC method to exact numerical results that we will ob- 
tain in a truncated Hilbert space.. To keep the system numerically affordable 
for the exact diagonalisation technique we consider a small 2x2 Mn02 lattice 
with periodic boundary conditions. The Mn sites are placed at integer coordi- 
nates {2i,2j) with i and j taking the values and 1 , while the oxygen atoms 
are at positions (2/ +l,2j) and (2/ , 2y + 1 ) . In order to limit to the maximum 
the dimension of the Hilbert space we restrict the degrees of freedom to those 
orbitals which are the most important for the physics of manganites. These 
are the eg 3d orbitals of Mn, namely the —y^ and 3z^ — orbitals, and the 
p oxygen orbitals which point toward Mn sites. For oxygens we restrict to 
for the {2i +l,2j) sites and py at the (2/, 2j +1) sites. These are the oxygen 
orbitals which bridge the Mn sites along the x and y directions. The system 
Hamiltonian is composed of several terms ; 

H^Hbare+Hhap+Hf}"+Hf"+H^ (19) 

namely Hhare which contains the one-particle energies of the orbitals, the 
hopping Hamiltonian Hhop which moves electron between neighboring sites, 
the Hubbard correlations H^", Hf" and for Manganese and Oxygen. 
The term is used because, applying exact diagonalisation, we truncate 
the Hilbert space by limiting the p orbitals occupation numbers between 1 
and 2. In the CC method, instead, we cannot truncate because this would 
destroy commutation relations. We have the possibility instead, in CCM, of 
choosing a high value of U in the Hubbard correlation , in conjunction 
with the oxygen part of Hi,are to effectively limit the oxygen p orbitals occu- 
pation numbers thus making the comparison, with the truncated model exact 
solution, possible. 

The bare Hamiltonian is 

Hhare = ^d.a '^gj.a,2i.2j'^gj,c,2i.2j + 

i^p - U„) ^ (/'L,,,^,,/,.„,2,+ 1.2, +/'L,2,.2;+l/'.v.a.2,.2;+l ) (20) 
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where the index takes the values — —y'^,3z^ — r^, with x^y be- 
ing in plane and z out of plane. The Mn one-particle energies are spin- 
dependent to take into account the mean-field exchange with the Mn t2g oc- 
cupied orbitals (xy,xz,yz)( whose degrees of freedom are discarded from the 
model). The oxygen orbitals term takes into account the Hubbard coeffi- 
cient —Up to compensate and favoring double and single occupations on 
oxygens. 

The hopping term is 

(21) 

where 

hz^-r2.,x^ hz^-r^.y = 1/2 
-fx2-y2,y = A2-v2,r = \/3/2 

The Coulomb intra-site repulsive interaction for Mn is made by a part 
for an electron pair on the same orbital, and another part for two different 
orbitals : 

^ .-^ grf.CT=+l 2,-.2j gj,a=-!,2i2j . ^ " 3z2-r2.(Ji,2,-,2j x2-y2,02,2ay 

(22) 

The Coulomb exchange for eg orbitals is 



Hf"=Jd y d\ ^ d\ ^ d, , d, , (23) 

. .'-I 3z2-r2,a22i2j x2-y2,oi.2,-.2; 3z2-r2,oi ,2,-.2; x2-v2.02.2l2; 



i.;,0i,02 

while the eg-t2g exchange is included as a mean-field term inside Hhare- 
Finally the oxygen Hubbard term is 

LUp [ n , n , +n , n , +2 

. . ' \ p,-.o=+2.2;+l,2; px.a=-2.2;"+l,2; pv.o=+2-2'-2;+I Py.o=-2.2i,2j+1 , 

(24) 

The contributions of the terms factored by Up, in the total Hamiltonian, 
is identically zero when we restrict the rip occupations between 1 and 2. 

To fix the free parameters of the model we use knowledge from our pre- 
vious work on manganites|91. Parameters are given in eV units. The effec- 
tive Slater integrals used in that work correspond, in the present model, to 
Ud = 6.88, U'l = 5.049, Jd — —0.917. The exchange with occupied polarized 
t2g orbitals gives a ~ 2eV splitting between e^(j^_i = 2 and £-j = 0, 
in the case of ferromagnetic alignement. We use a hopping f = 1.8 taken 
from our previous work|9 |. The parameter Ep controls the amount of charge 
back-donation from oxygen to manganese. The predominant O 2p character 
of doped holes found in manganites ifTOl corresponds to a value which 
raises the bare oxygen orbitals energies above the bare Mn ones. The value 
of Ep influences the average occupation of the orbitals. These occupancies 
match the ones found in the previous works for a value Ep ~ 2. 
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4 Discussion 



To find the CCM ground state and determine the resolvent equation we have 
adapted our Hilbert + + ifTTl |9] code. This code was originally created to 
calculate x-ray spectroscopies of small strongly correlated clusters by ex- 
act diagonalisation. It implements a second quantisation representation of 
operators and determinants. We have implemented automatic computing of 
commutators and automatic extension of the excitations set for CCM and for 
our sec method. 

The exact diagonalisation and the Lancsoz tridiagonalisation for spectra 
calculation are performed with Hilbert + +. The code generates the Hilbert 
space by applying several times the Hamiltonian on a vector basis which 
is beforehand initialized with a seed state. In this seed state, the occupied 
spinorbitals are all the oxygen ones and, for ferromagnetic alignement on 
the Mn sites, all the spinorbitals — with spin o = +1/2 . This state is 
named, in the rest of this paper, nominal reference configuration. The con- 
figurations having one or more oxygen sites unoccupied , are discarded in 
the exact calculation. With this limitation on the configurations, the gener- 
ated Hilbert space growths up to a dimension which is slightly less than 7 
millions. 

To reproduce with the CCM method the exact calculation done on the 
truncated space we set Up as high as IQ^eV . 

We show in figure [T] the convergence of the CCM energy as a function 
of the number of symmetrized eccitations contained in the S operator when 
we take the nominal reference configuration as reference state The CCM 
energy converges, for the nominal reference, to the first excited eigenenergy, 
given by exact diagonalisation, above the ground state. We have analysed the 
ground and the first excited states that we obtain by exact diagonalisation. 
The largest component of the first excited state is found to be the nominal 
reference state. This explains why the CCM method, which takes this state 
as reference, converges to this eigen-state. The ground state, instead, has a 
different symmetry. We find that there are four components which have the 
largest factor and each of this component is obtained rotating the e^ electron, 
on one of the four Mn sites, from the — r^ orbital to the —y^ one . More 
in details, the ground state has the same symmetry of the state |<l>o) given by 

= T{-'^Y'"^d\ , d, , lOo) (25) 

I / ' x2-,.2.o=l/2.2,-.2; 3z2-,2,o=,/2.2i.2; I ^ ' 

l-J 

This state cannot be obtained starting from the nominal reference with the 
CC method because it has completely different symmetry properties. Notice 
for example that a 7t/2 rotation around the center of the cluster gives a factor 
1 if applied on the nominal reference but the same rotation gives a factor — 1 
when applied on j't'o)- 

To force the CC method to converge to such a symmetry state, a possible 
solution would be using the multi-reference CC method. We have not imple- 
mented this method, which requires an a-priori knowledge of the solution, 
because our work is focussed on our Spectral Cluster Method that, as we will 
see in this section, is able to detect these states exploring selected spectral re- 
gions. To test further the capabilities of CCM to converge to the true ground 
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state we have allowed a possible convergence to |<I>o) symmetry by using a 
reference state |<l>o) of lower symmetry : 

|*o) , d , , lOo) (26) 

I ' .v2-y2 .0=1/2,0.0 3z2-,-2 .0=1/2.0.0 ' 

We show in figure|2]the convergence of the CCM energies for Up — 5eV 
using two different choices of the reference state : the high symmetry state 
|<t>o) and the lower symmetry j'l'o) ■ The CCM solution for reference state 
|5>o) converges to a lower energy than the one obtained with the nominal 
reference state. We cannot compare this calculation, done with Up = 5eV , 
with the results of exact diagonalisation because the small value of Up gives 
access to a larger Hilbert space which is computationally more expensive. 
On the other hand for an high value of Up — lO^eV, when comparaison with 
exact diagonalisation is possible, we have not been able to obtain the ground 
state starting from the low symmetry reference state |5>o)- We think that 
this difficulty can be explained in the following way : the lower energy of the 
|<I>o) symmetry state is due, in the CCM equations, to a kind of bridges, made 
of operators which link the components of |<I>o) to each other. These bridges 
are created when an excitation operator, which composes S, is transformed 
by the Hausdorff commutation expansion into another excitation operator 
which will subsequently enter S, and so on. For the particular symmetry of 
|<I>o) to be obtained from j'l'o), these bridges must be long enough to trans- 
form one component into another The problem of using a high value of Up 
is that for every pair of excitation operators which both create a hole on the 
same oxygen site, a new term coming from their product, will appear in the 
residue containing two holes that site. This will necessitate a new higher 
order excitation to be subsequently included in S whose contributions will 
cancel the product of the two operators. This because the very high value of 
Up forbids double hole occupancies on oxygen sites. The need of account- 
ing more operators requires more iterations. During these iteration the |<I>o) 
symmetry is non favorable and our procedure converges to higher eigenval- 
ues. The CCM wavefunction corresponding to the true ground state becomes 
energetically favorable for a number of eccitation operators about 40. The 
use of a multireference could have been used to force a particular symmetry. 
This procedure would have been feasible for the small system that we have 
treated in this work, because we can know the ground state symmetry from 
the exact solution. For larger systems, however, even if one could a-priori 
know the correct symmetry, the number of determinants in the multirefer- 
ence state grows exponentially with the size of the system. Moreover, the 
Newton method, used for solving the CCM equations, does not guarantee 
that the lowest energy solution will be found, because this method allows to 
follow just one solution which might be not the good one. 

The spectral method that we present in this work allows instead to ex- 
plore, focussing on selected spectral regions, a larger set of solutions which 
are observed as resonances. 

We show in figure 3 the spectra for the first excited eigenvalue at = 
2ey, Up = 10^ eV, considering a probe operator D^f which induces a crystal- 
field rotation in the Cg space, on one Mn site : 

Dcf = d\ ^ d ^ , (27) 

■' .v2-v2 .0=1/2.0.0 3z2-r2.o=l/2,0,0 
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The SSCM equations reproduces well the exact diagonalisation results. The 
spectra shows a peak at negative energy. This is the ground state which was 
not accessible starting from the nominal reference state but it is visible as 
a resonance in the SCCM equations. The SCCM residues, to expand the R 
operator, have been calculated fixing cOr at zero because crystal field excita- 
tions are found at low energies. The SCCM spectra has been calculated with 

Figure 4 shows the spectra for the same initial state, but for a probe oper- 
ator which transfers charge from a oxygen site to a neighbouring manganese: 

Dcf = d\ , p . (28) 

The SCCM spectra has been calculated considering two energy windows : 
one around the charge transfer peak using co, = ?>eV,N'' = 8 x 10^ , and an- 
other window around the ground state, using co^ = — 0.5ey,A'^'" = 10*. 

Figure 5 show the same spectra calculated for a non-truncated Hilbert 
space, using Up = 5eV. Comparison to exact calculation is not possible 
in this case but we can see that the most important features are preserved, 
namely the charge transfer peak, the ground state peak and the peak due to 
the overlap with the initial state at zero absorbed energy. The convergence of 
the spectra in this case of low Up is easier and the spectra can be calculated 
with only one energy window, using co^ = TieV . The graph shows two curves, 
one calculated with A^'' = 10^, where the charge transfer peak is already in 
place, and another done at a higher value of AT'' = 7 x 10^ which is neces- 
sary to have a proper convergence on the ground state peak at ~ —O.leV. 
The different behaviour for the two peaks can be seen as a consequence of 
the non-locality of the ground state derived from |4>o) symmetry. The non- 
locality impUes a larger set of terms entering the resolvent simi. 



5 Conclusions 

We have appUed CCM equations to a strongly correlated lattice in the case of 
strong departure from the reference state. We have developped the spectral 
coupled cluster equations, by finding an approximation to the resolvent oper- 
ator, that gives the spectral response for the class of probes that are writable 
as products of creation operators 

We have applied the method to a Mn02 plane model for a parameters 
choice which makes the ground state particularly difficult to find with the 
CCM equations because of its peculiar symmetry which corresponds to a 
non-nominal reference state. We have shown that this state can be spectrally 
observed using SCCM equations by probing a CCM solution for the nominal 
reference state. In this case one observes a negative energy solution which 
corresponds to the true ground state. We think that CCM and SCCM equa- 
tions have a strong potential, for strongly correlated lattices, not only for the 
study of the ground state but also for all those excitations that can be repre- 
sented by a resolvent operator R that can be written as the sum of locahsed 
terms. 
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Figure 1 : CCM versus exact diagonalisation of the truncated Hilbert space space. 
For the CCM method we apply an effective truncation by setting Up as high as 
10^ eV. The convergence of the CCM energy is shown as a function of the number 
of symmetrized eccitations contained in the S operator. The nominal reference 
generates the first excited state while the ground state has the symmetry of |4>o)> 
equation [2T 
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Figure 2: Convergence of the CCM energies for a non-truncated Hilbert space 
(Up ~ 5eV). The nominal reference generates the first excited state while the lower 
symmetry reference state |<l>o) gives the ground state. 
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Figure 3: Spectra for the out-of-plane to in-plane probe Dcf applied on the first 
excited eigenvalue, Up = 10^ eV (truncated hilbert space) 
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Figure 4: Spectra for the o charge transfer probe Dch applied on the first excited 
eigenvalue of the truncated problem {Up = 10^ eV ) 
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Figure 5: Spectra for the charge transfer probe Dch applied on the first excited 
eigenvalue of the non truncated problem {Up = 5eV ). 
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